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We study the formation and evolution of an interconnected string network in large-scale field- 
theory numerical simulations, both in flat spacetime and in expanding universe. The network 
consists of gauge (7(1) strings of two different kinds and their bound states, arising due to an 
attractive interaction potential. We find that the network shows no tendency to "freeze" and 
appears to approach a scaling regime, with all characteristic lengths growing linearly with time. 
Bound strings constitute only a small fraction of the total string length in the network. 
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. It has been recently realized that fundamental strings and other string-like objects, such as 15-strings, can have 
' astronomical dimensions and play the role of cosmic strings.^ Observing these objects in the sky would provide the 
D , most direct test of superstring theory. Both types of string are naturally formed in the course of brane-antibrane 
pLn ' annihilation at the end of brane inflation 0, [3, [j, d, 0] ■ Fundamental (F) and £)-strings produced in the aftermath 
of this annihilation can form {p, q) bound states combining p F-strings and q D-strings. As a result the strings are 
expected to form an interconnected FZ)-network 0, [6|, with different types of string joined in 3- way y-type junctions. 

Similar string networks can also be formed in field theory; a simple example has been recently given by Saflin 0]. 
His model includes two Abelian Higgs models, with an additional coupling between the Higgses. The model has a 
broken U{1)a x U{1)b symmetry, resulting in two types of string, and the coupling is chosen so that A and _B-type 
strings are attracted to each other and can form (p, q) bound states. An even simpler example is the usual Abelian 
Higgs model. With a suitable choice of the Higgs and gauge couplings, corresponding to the type-I regime, this model 
, allows stable strings with arbitrary windings, which can be joined in 3-way junctions Strings with extreme type 

I properties can also be formed in models with SUSY flat directions 
] String networks with 3-way vertices can also arise in models with symmetry breaking of the kind G — > [lol | , as well 
^ . as in non- Abelian field theories with several types of strings corresponding to non-commuting symmetry generators 
' Ii3- latter case, when two non-commuting strings cross, a third string starts stretching between them, 

resulting in two y-j unctions. 

The evolution of cosmic string networks has been a subject of much recent discussion and debate.^ Early work on 
Za-strings using a simple analytic model, suggested that at late times the characteristic scale of the network, 
C = it^/PsY^'^i exhibits scaling behaviour, 
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Ps/p ^ G^ih\ (2) 



Here, ps and p are the average energy densities of the network and of the universe, respectively, /i is the string tension, 
and 7 is a constant . A similar model was later used for networks having N strings joined at each vertex p^ . 

The magnitude of 7 in Eq. ([T]) depends on the rate of energy loss by the network. In models, the vertices can 
carry an unconfined magnetic charge. The energy of the network is then efficiently dissipated by gauge field radiation 
from these magnetic monopoles. Another energy loss mechanism is the formation of closed loops and of small nets 
disconnected from the main network. In the absence of magnetic charges, and if loop and net formation turn out to 
be inefficient, the remaining energy loss channel is the gravitational radiation. In this case, the analysis of [15] gives 
7 ~ G/i, and Eq. ([2]) gives Ps/ p ^ 1/G/i ^ 1, indicating that the string network becomes so dense that it dominates 
the universe. 



^ Higher-dimensional D-branes with all but one dimension compactified will also appear as stringlike objects from a maeroscopic point of 
view. 

^ One of the key questions is whether or not the network gets entangled and "freezes", in which case it would eventually dominate the 
energy density of the universe [Tsl. Il4l|. However, most recent work points to scaling evolution. 
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More sophisticated analytic models have recently been developed, aimed directly to describe cosmic superstrings. 
These models allow for several types of string with different tensions and use the velocity-dependent one-scale model 
of string evolution [13, [3 (see also [l^). The models make somewhat different assumptions about the physics of F- 
and I?-string interaction. Tye, Wasserman and Wyman flT'l assume that when F and D strings meet and "zip" to 
form a bound _FZ?-string, the excess energy is released in the form of high-energy particles. If this picture is correct, it 
provides an important additional mechanism of energy loss by the network. They also assume that the entire network 
is characterized by a single length scale £,{t). Avgoustidis and Shellard fisl assume, on the other hand, that the energy 
released in the zipping process goes to increase the kinetic energy of strings, and thus remains in the network, and 
allow different length scales for different string types. Assuming that the rate of energy loss to loops from networks is 
about the same as that from "ordinary" strings, both models predict scaling evolution, with 7 ~ 1, with energy about 
equally distributed between (1,0), (0,1), and (1,±1) strings, and with negligible amount of energy is higher-(p,q) 
strings. 

Network evolution has also been studied in field-theory numerical simulations.'^ Spergel and Pen |20] and later 
Copeland and SafSn used a non-linear sigma-model to simulate non-commuting string networks. Hindmarsh and 
Saffin [HI performed a full field theory simulation of global strings. In all this work, the network was found to scale 
with 7 ^ 0.1 — 1, indicating efficient damping. The models used in [l6l. [20l. [2I1] have some important differences from 
superstring networks. First, all types of string in these models have the same tension, while in an FD-network the 
tensions of all (p, q) strings are generally different. Second, the global symmetry breaking models used in [l6l. [20I [2H 
allowed for an additional energy loss mechanism - the radiation of massless Goldstone bosons - which is known to be 
rather efficient. On the other hand, superstring networks are expected to have only gravitational-strength couplings 
to massless (or light) bosons.^ Hence, there is a danger of string domination, unless the network can efficiently lose 
energy by loop or small net production. 

A field theory simulation of an interconnected string network has been recently developed by Rajantie, Sakellariadou 
and Stoica [1^ . They used a model of interacting scalar and gauge fields similar to the Saffin's model |7i] , which allows 
two types of string and a spectrum of bound sates. The dynamic range of their simulations was not sufficient to reach 
any conclusions about the statistical properties of the network and its scaling behaviour (or lack thereof). The main 
focus of Rajantie et. al. paper is on the effect of the long-range interaction induced by the Goldstone field in models 
where one of the two U{1) symmetries of the model is global. They find that this interaction disrupts the string bound 
states in the network. This result is probably of little relevance for superstring FZ)-networks, since, as we already 
noted, superstring interactions are expected to have gravitational strength and thus have little effect on network 
dynamics. 

In this paper, we have developed a new network simulation, which we believe to be closer to a "realistic" superstring 
network. We used a field theory model with two types of gauge strings and adopted Saffin's [3] interaction potential 
to ensure that the strings form bound states. The details of the model and of the simulation are given in the next 
section. The results are presented in Section 3. Our conclusions are summarized and discussed in Section 4. 

II. SIMULATION DETAILS 
A. The model 

Safhn's model of interacting strings^ is defined by the Lagrangian [3] 

c = \D,^\' + \v,iA^ \fi^ \ti ~ vm, m. (3) 

Here, cf) and tp are complex scalar fields, charged with respect to and gauge fields, respectively, 

F>f, = dfj,- leAfj,, 'D^ = d^- igBf,, (4) 

Ff,u = d^,A^ - J^^^ = d^,B^ - d^B^,, (5) 



^ Earlier simulations, using a simple model of straight strings joined at vertices, were performed in 15] for a Z3 network and in Il2ll for 
non-Abelian strings. 

* The possibility of cosmic superstrings having stronger than gravitational couplings to massless Ramond-Ramond fields has been recently 

discussed by Firouzjahi [23 |. 
^ Similar models have been studied in relation to composite defects l24l I25I . |26| . l27tl . 
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v{\<p\M)^^ 



v\f + ^m'~ vir - -(101^ - vDm' - vD- 



(6) 



Without the last term in the potential, the model describes independent A- and S-strings. Bound states are formed 
if the parameter k is chosen in the range 0] 



< K < ^^/XaXi 



(7) 



In this paper we shall not attempt to explore the full parameter space of the model and consider only the special 
case where the strings are in the Bogomol'nyi limit, 

Xa = 2e^ As = 2.g2. (8) 

We shall also set e — g and tja — rfs- With standard rescalings, the parameters of the model can then be reduced to 

VA=VB = 1, (9) 



and the range of k in Eq. ([7]) becomes 



^Xa = Ixb 



< K < 1. 



1, 



(10) 



(11) 



In most of our simulations we used the value k = 0.9. Table U gives the corresponding string tensions, as well as the 
binding energies (per unit length of string), which are relatively large. We have also included in Table [J the values 
for a larger k = 0.95. 



(m, n) 


K = 0.90 
fi A/i 


K = 0.95 
/i A/i 


(1,0) 
(1,1) 
(2,1) 


0.793 - 
1.278 0.308 
1.798 0.581 


0.728 - 
1.133 0.323 
1.560 0.624 



TABLE I: String tensions (^) and bounding energies (A/i) for (m, n) type strings, calculated for k = 0.90 'T] and k = 0.95. 



B. Numerical setup 



Our aim was to perform real-time lattice simulations of model ([3]), for as long a time as our facilities allowed us. 
In order to do this, we obtained the equations of motions from Q, 
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(12) 



where we have made a gauge choice Aq ^ Bq = and assumed a flat FRW spacetime written in conformal time, 

ds^ ^a2(T)(dT2 -dx^). (13) 

Overdots in Eqs. (fT2|) stand for derivatives with respect to r. 

There is a well-known problem in such simulations: the string core has a fixed physical width, whereas the distance 
between lattice-points grows with the expansion (a oc r in radiation era, and a oc in matter era). As a result the 
string width quickly drops below the resolution threshold of the simulation. Here, we have adopted the approach used 



4 



in [28|, [291, [SO, l31| , in which the equations of motion are artificially modified to have the string width growing with 
the expansion, so as to be able to simulate them throughout the evolution. 
Following [Sil] , the equations of motion ([T2|) can be written as: 



a . 



Foj - 2(1 - s)-diF,j = -2a^"e lm[(f)* Dj(j)\ (14) 

and likewise for ij} and B^. Here, s is a parameter that controls the growth of the string width, with s = 1 being the 
true value. As earlier work has shown [30l . [sij . there is little difference in string dynamics for different values of s. 
For the remainder of this work, we shall set s = (that is, the string has constant comoving width). 

We discretized the modified equations of motion ([T^ on a lattice using the standard lattice link variable approach 
[3^ and performed the simulations on the UK National Cosmology Supercomputer [ssj . The simulation box consisted 
of 512'^ lattice points, with periodic boundary conditions. We chose Aa; = 1.0 and Ar — 0.2 trying to maximize the 
dynamical range of the simulation. 

Scaling evolution regime is expected to be an attractor, and indeed earlier work has shown jsi], [s^, [H, that this 
regime is approached from a wide range of initial configurations. Nonetheless, constructing initial conditions for this 
kind of simulation is a nontrivial task. The challenge is to find some initial configuration that leads to scaling as fast 
as possible, in order to maximize the dynamical range. 

For the results presented here, we used the following procedure: we set all gauge fields and gauge field momenta to 
zero, also set the scalar field momenta to zero. The magnitudes of the scalar fields is set to rjA and rjB respectively 
(in the present work rjA — i]b — 1), and their phases are chosen randomly. This configuration was then smoothed 
out by averaging over nearest neighbours, and this procedure repeated 20 times, in order to get rid of some excess 
initial gradient energy. Note that this initial configuration satisfies the (lattice) Gauss law, and due to the lattice-link 
variable procedure, the Gauss law is guaranteed to be satisfied throughout the simulation. 

The initial configuration described above corresponds to a highly excited state, so in order to remove the extra 
energy from the system, a fake constant damping was applied, 7 = f = 0-2, until time r = 32. From then on, we 
performed simulations for flat spacetime (7 = 0), radiation (7 = 1/r) and matter (7 = 2/t) eras. Only times r > 64 
were used for this analysis. 

In order to automatically detect the strings in the simulation and compute their length, we calculated the net 
winding of the phases around plaquettes. One can then trace the string following the winding and estimate the 
lengths. The length L of a string that crosses n plaquettes is estimated as L — n ■ Ax. 

An AB string can be traced by the sites where both A and B phases wind. Unfortunately, there are two drawbacks 
of this procedure: on the one hand, there are accidental crossings of A and B strings, i.e., lattice sites where an 
A string and a B string simply cross and follow their way, without any attempt to form a bound state. On the 
other hand, there are places where inside a clear segment of AB string, A and B phases do not wind in exactly the 
same plaquettes, but there is a slight displacement (we will show an example in the following Section). In order 
to circumvent these problems, we have excluded accidental crossing by setting a minimum size of AB string {Imfn)^ 
and also excluded gaps in AB segments (due to an occasional "displacement") shorter than a given distance {d^^^). 
This means that in the process of detecting AB segments, all segments of size I < will not be considered as 
legitimate segments. Also, if the separation between two segments is d < d1^^, those two segments are considered as 
one. The total length of AB string is not considerably affected by this process; the main difference corresponds to a 
more realistic count of the number N of ^-junctions formed, that is, it helps in not overcounting AB-segments and 
y-junctions. There is no obvious way of determining the values of /j^^ and d^^^, and we have chosen them to be 
^mfn = 3 and d^^^ = 5 by trial and error, and inspection of the results. With this choice, the corrected value of N is 
roughly a factor of 4 smaller than one would get from the raw data. 

As we already mentioned, our choice of parameters was largely motivated by the effort to increase the dynamic 
range of the simulation. For example, we allocated only a few lattice points per string thickness. As a result, our 
discretized representation of the field theory string solutions is not particularly accurate. For a rapidly moving string, 
this may result in spurious damping, with the kinetic energy of the string being dissipated into particles [IJ- Moore 
et. al. [201 performed numerical tests to determine the optimal choice of the lattice spacing Aa: and concluded that 
Ax = 0.5 is close to the maximal value that still accurately represents the string dynamics. For larger values of Aa;, 
they observed a significant increase in spurious particle emission from oscillating strings. 

In the present paper, our focus is not so much on the dynamics of oscillating loops, as it is on the overall charac- 
teristics of the large network. In fact, in order to observe oscillating loops in field theory simulations Aa; should be 
rather small (even smaller than 0.5), leaving us with a tiny dynamical range where large network properties would 
be impossible to study. We therefore pushed the parameter values a bit further and used Aa; = 1.0, with the hope 
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that these properties will not be strongly affected. We performed several tests by simulating the system with different 
Aa; (0.5 < Ax < 1.0), smoothness of the initial configuration (smoothing between and 30 times), the value of the 
initial fake damping (between and 0.5) and the length of time in which the fake damping was active {0 < t < 50). 
The qualitative results for all characteristic lengths of the network were the same in all cases, with the actual values 
agreeing within 10 — 20%. (The only exception is the length of A and B segments, as discussed in Sections [ill CI and 

nmi) 




FIG. 1: Picture of a typical simulation of a (p,q) network. The green and blue colours correspond to A and B strings respectively, 
whereas the red colour shows the AB segments. It can be clearly seen how Y junctions are formed all over the simulation. 



III. RESULTS 

A. The network correlation scale ^ 

Our simulations were performed in flat spacetime and in the radiation and matter eras. In all cases an interconnected 
network was formed with A-, B- and (1,1) ^i?-strings. No higher-(p, g) strings were observed. A representative 



6 



snapshot of the network is shown in Fig. [T] Throughout the evolution, the network is dominated by one large 
("infinite") interconnected net, comprising more than 90% of the total string length, as seen in Fig. [21 




FIG. 2: Picture of a typical simulation of a (p,q) network. This picture represents the same configuration, but we show how 
most of the string length (typically more than 90%) forms an interconnected network (A and B strings in the main network 
shown in green and AB segments in red) and there are only a few loops that do not belong to the main network (blue). 

As mentioned in the previous section, within a single AB bound state, A and B strings can be displaced by a 
single lattice point, making the code decide that it is in fact two separate segments. Figure [3] shows a fragment of 
the simulation box, with a somewhat long AB string depicted. Those accidental displacements should not be taken 
into account, and with the help of the parameter d^^^ (introduced earlier), the displacements are reassessed, and the 
segment is counted as one. 

The overall length scale of the network can be defined as usual, 

^={V/Ly/^ (15) 
where V is the volume of the simulation box, L is the total length of string. 



L — La + Lb + Lab, 



(16) 
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FIG. 3: Fragment of the simulation box showed in Fig. [T] Green and blue correspond to A and B string, whereas the red 
colour corresponds to an AB segment. A long AB segment can be seen in the picture, but at some points the A and B string 
miss each other by just a lattice point. With the help of the parameter d^^^ those accidental displacements are accounted for, 
and long segments such as the one in the figure are counted as one. 



and La, Lb and Lab are the lengths in A, B and AB strings, respectively. ^ gives the typical distance between 
strings in the network. 

Throughout this paper we shall use comoving length scales. The corresponding physical lengths, which will be 
denoted by superscript (ph), can be obtained by multiplying with the scale factor a(r), e.g., 

e(^"')(r)=a(r)C(r). (17) 

For our models, the scale factor has the form 

«W = (r/ro)", (18) 

with To = const and n = 0, 1, 2 for flat spacetime, radiation and matter eras, respectively. The physical time t and 
the horizon distance £h are given by 



£h = a{T)T = {n + l)t. (20) 

The simulation results for ^(r) are shown in Fig. [31 where each graph represents an average over 10 simulations. 
Remarkably, the graphs for the flat, radiation and matter regimes are almost identical. They show a nearly linear 
dependence, 

e(r)=ar + C(«), (21) 

with a « 0.15 (see Table lll| . Toward the end of the simulation, the two terms in Eq. (PT|) are comparable to one 
another. If the linear dependence extends to much larger values of r, the constant term eventually becomes negligible, 

f « ar, (22) 

and the corresponding physical length grows proportionally to the horizon, 



^'^P''\t) « {n + l)at = a£H{t). 



(23) 
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FIG. 4: The correlation length ^, as defined in Eqn. 1151 averaged over 10 simulations for flat (dashed black), radiation 
(continuous red) and matter (dotted blue) regimes. Note that all three cases exhibit an approximately linear behaviour, nearly 
independent of the regime. 





Flat 


Radiation 


Matter 


«(C) 


0.14 


0.15 


0.15 


a A Ha) 


0.21 


0.22 


0.21 


cy-AB (Iab) 


0.03 


0.07 


0.08 


Q-N Hn) 


0.21 


0.28 


0.30 



TABLE II: The values of the linear growth coefficients of different lengths, as defined in the body of the paper, obtained by 
fitting the simulation data. Note that with our parameter values a a = Q-b- 



B. Bound strings 

An important characteristic of the network is the fraction of total string length in bound {AB) strings, 

fAB = , ^Y\t • (24) 

La + Lb + Lab 

The simulation results for Jab are shown in Fig. [5l We see that Jab remains nearly constant, at the value 
0.01 < Jab < 0.02 in radiation and matter eras, and somewhat lower for flat spacetime. Hence, bound strings 
constitute less than 2% of the network. This is in conflict with analytic models (TtI. Il8l| predicting that the energy of 
the network should be more or less equally divided between A, B and AB-strings. The main source of the discrepancy 
is the assumption made in [13, [3] that crossings of A and B strings typically lead to the formation of relatively stable 
bound AB segments of length ~ £,{t). Visual inspection of the simulation movies suggests, on the contrary, that 
formation of bound segments by intersecting A and B st ring s occurs rather infrequently, probably when the relative 
velocity of the colliding strings is sufficiently small (38|, [39l |. Even when they are formed, the Ai?-segments easily 
"unzip" as the free A and B ends pull in different directions and do not usually last for more than a Hubble time. 

Even though the fraction of bound string is small, the interaction of A and B-strings has a significant effect on 
the network evolution. To quantify this effect, we ran some simulations with the same initial conditions for A-string 
fields as before, but with i?-string fields starting in their ground state, so that no B strings are formed. A-strings 
then evolve as ordinary U{1) strings, and their characteristic length 

U = {V/LaY^^ (25) 

exhibits a linear dependence on r with ua ~ 0.29, as seen in Fig. [6l (This value is in agreement with earlier U{1) 
simulations by Vincent et. al. [13] and Moore et. al. [33|.) The same quantity calculated with B-strings present gives 
ua ~ 0.22, so the growth of S,a is significantly slower than it would be if the two kinds of string evolved independently. 
The average (comoving) length of AB-segments is 



Iab = Lab/N, 



(26) 



9 



0.1 



CD 
< 



0.01 



Flat 

Rad 

Mat 



0.001 

80 100 120 140 160 180 200 220 240 260 
time 

FIG. 5: Fraction of total string length in bound strings, for flat, radiation and matter regimes. The percentage is fairly low; 
between 0.5-1% for flat, and 1-2% for radiation and matter cases. 
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FIG. 6: Comparison of (only in radiation era, for clarity) between simulations points and without B-strings. 



where N is the number of AiJ-segments (2A^ is the number of F-junctions where the three types of string meet). 
Fig. [7] shows that the evolution of Ij^^g is approximately linear, 



Iab ~ CXABT + I'AB 



(27) 



The coefhcient aAs (see Table is nearly the same in radiation and matter eras (within 10%), 0^4^'''""'*'' w 0.07, 

and is significantly smaller in flat spacetime, a^^"*^ w 0.025. The shorter bound segments in flat spacetime are 
probably due to larger string velocities. 

The typical distance between y-junctions is given by 



Cat - iV/2N)^/\ 
Once again, we find approximately linear evolution (see Fig. [5]), 



(28) 



(29) 



with aN ~ 0.3. 
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FIG. 7: The average length of bound segments Iab, as defined in Eg. 1261 for flat, radiation and matter regimes. The evolution 
is approximately linear, with slope given in Table [III 
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FIG. 8: The average distance between y-junctions, as defined in Eq. 1281 for flat, radiation and matter regimes. The 
evolution is approximately linear, with slope given in Table HIl 



C. Non-scaling of A and B segments 

All results presented so far are consistent with scaling evolution, with all characteristic length scales of the network 
growing proportionally to the horizon. However, the average comoving lengths of A and B segments, 

Ia = La/N, Ib = Lb/N, (30) 

do not exhibit scaling behaviour. In fact, Fig.[5]shows that, rather surprisingly, these lengths approach nearly constant 
values, 

Ia ~ Ib ~ const. (31) 

The corresponding physical lengths grow proportionally to the scale factor. 

Throughout the simulation, the segment length I a remains much greater than the correlation length ^a- The 
correlation length grows with time, but is still about an order of magnitude smaller than Ia at the end of the 
simulation (and similarly for Ib and ^b)- This suggests that A (and B) segments have the shape of random walks of 
step ^ ^A and end-to-end distance ~ ^n- The length of the segments is then 



(32) 
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FIG. 9: The average length of A-segments Ia is shown for flat, radiation and matter regimes. In all three cases, Ia is nearly 
constant at late times. The behaviour of Ib is essentially the same. 



Physical time 


Damping (7) 


At 


< t < 32 
32 < t < 272 
t > 272 


0.5 
2.0 

2.0/t (rad.) 


0.2 
0.5 
0.2 



TABLE III: Values of the parameters in simulation with a second damping period, introduced to achieve scaling of the segment 
lengths I A and Ib- 

We have verified that the ratio £,%/£,aIa = (.a/'^^n is indeed approximately a constant of order 0.3, within 10 — 20%. 

If the length of the segments I a continues to grow slower than their end-to-end distance £,n , the two lengths will 
eventually become comparable to one another, with segments becoming more or less straight. This could mark 
the beginning of the true scaling regime, where all the characteristic lengths of the network have the same order of 
magnitude and grow proportionally to r. In any case, the evolution laws ()2ip .(|27p. (j31|) cannot continue indefinitely and 
must stop at some r = r, . The situation here may be somewhat similar to that with "ordinary" (not interconnected) 
strings, where scaling of the characteristic length ^{t) is quickly established, but scaling of the small wiggles on long 
strings and of closed loops is reached only after a long transient period [41,, ,4^, [4^, ,4^] . We shall refer to the evolution 
at T < T» as the transient scaling regime. 

D. Towards a true scaling regime 

We attempted to shorten this transient regime, or avoid it altogether, by increasing the duration of the initial 
damped period. This has the effect of increasing N and decreasing the ratio Ia/£,n ^ (Ca'/Ca)^ in the initial state 
(right after damping is turned off). Getting to the end of the transient regime by the end of the simulation proved to 
be a difficult task, for which we had to push the limits of the stability of the code by using a rather large time step, a 
rather large damping coefficient, and evolving the system beyond the half-light-crossing time of the simulation box. 

We performed a set of 20 simulations with the values of time step At and the damping parameter 7 shown in 
Table IIIII The difference from the rest of our simulations is in the second period of damping imposed after the first 
period (which is common to all simulations and is used to relax the system from the highly excited initial state and 
allow string formation). The extra damping period has a time-discretization of At = 0.5 (keeping Ax = 1), and the 
numerical stability is achieved because the damping is also rather high (7 = 2.0). 

Figure [TUl shows the length measures ^a, £.n, Ia averaged over the 20 simulations. All three lines seem to approach 
a linear behaviour. The slopes of ^a and agree with the ones obtained earlier, without an extended damping 
period. Other quantities characterizing the system, such as oab, Jab and the percentage of string length in loops, 
also agreed with the ones for the simulations without extra damping. The fairly linear behaviour of Ia is the main 
change. A linear fit to the evolution of I a gives 




I A ~ aLAT + Ia , 



(33) 
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with aLA ~ 0.8. Extrapolating this behaviour to large r, we expect a scaling regime in which Ia is a few times larger 
than ^ and ^jy. However, because of our extreme choice of parameters and large error bars in Z^, and because the 
change in Ia during the linear regime is relatively small, this conclusion should be regarded as tentative. 




FIG. 10: The length measures ^a, Cat, Ia averaged over 20 simulations performed with an extended damping period, as explained 
in Section HIT CI The thick lines corresponds to the average, whereas the shaded regions correspond to l-cr statistical errors 
over the simulations. 



E. Loops and small nets 

We observed the formation of small independent nets a few times in the course of the simulation, but these occasions 
were rather rare, so AB-strings belonged predominantly to the infinite network, with at most one small net in addition. 
The fraction of the total string length in disconnected A- or i?-loops is shown in Fig. [TlJ We see that this fraction 
remains nearly constant, at a value ^ 0.03 — 0.05 in radiation and matter eras and /l ^ 0.08 in flat spacetime. 
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FIG. 11: Percentage of string length for both A and B strings in loops that do not belong to the main network. 



The length distribution of independent loops is plotted in Fig. 1121 averaged over 20 simulations in radiation era. 
The figure uses logarithmic binning and shows the distribution of loops at four different times. We see that some of 
the loops are fairly large, with length much greater than the inter-string separation ^. Such loops should arguably 
be regarded as part of the infinite network. Indeed, if the diameter of a loop is larger than a few times £,a (or ^b) 
for a given time, the loop is very likely to reconnect to the network. If we exclude loops longer than 6^a (J^^a), 
the percentage of string length in loops drops to 1 — 2% (0.5%). These values are in agreement with previous field 
theoretical simulations where a loops were found to account for a few percent of the total string length [dO] ■ 

Simulation movies show that loops that decouple from the network do not oscillate as they would in Nambu-Goto 
simulations, but rather shrink and disintegrate. This could be expected, since it is well known that in order to observe 
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FIG. 12: (Left) Distribution of independent loops for the radiation era, averaged over 20 simulations, for four different time 
steps. The figure shows that there are rather large independent loops formed during the simulation, and that the peak of the 
distribution tends to move towards larger sizes with time. (Right) Loop distribution averaged over 20 simulations in radiation 
era. The distribution appears to follow Eqn. [34] rather well, suggesting that it is scaling. 



oscillating loops in field theory simulations one would need very large loops and very small values of the lattice spacing. 
Additional damping, through particle emission from loops, may be due to the presence of short-wavelength string 
excitations in the initial conditions, as indicated by the numerical results of Refs. [l^H^. The short lifetimes of the 
loops in our simulations explain, at least in part, the relatively small amount of string in loops, as compared to the 
Nambu-Goto simulations. 

In any case, we checked the loop distribution in our simulations for scaling behavior. On dimensional grounds, the 
number density of loops per unit length interval in a scaling network should have the form 

'^(t)^t-'f{l/t). (34) 

Defining /S.Ni as the number of loops between sizes U and U+i and AZ^ = li^i ~ h, we plot in Fig. [T^lthe (binned) 
quantity ^^f^, obtained for 5 different times and averaged over 20 simulations. In a situation where the loop 
distribution scales, all the lines should line on top of one another. We see that the graph in Fig. [T^ does indeed 
exhibit scaling. 

F. Effect of higher binding energy 

We ran some simulations with a greater value of the string binding parameter, k = 0.95, in order to try to understand 
how sensitive our results were to the precise value of k. The main difference encountered in these simulations was 
the fraction of bound string Jab increased by a factor ^ 1.5 (see Fig. [T5)) . Our analysis showed that the effect on 
the correlation lengths ^ and on the loop distribution was very small. Actually, the values of the coefficient of the 
different length measures where the same as in the k — 0.9 case, the only difference being a smaller for this case 
(0.16,0.22,0.25 for flat, radiation and matter eras, respectively). 

IV. CONCLUSIONS AND DISCUSSION 

We performed numerical simulations of the formation and evolution of string networks using the Saffin's model 
for interacting strings. The model has two types of Abelian gauge strings, A and B, with an attractive interaction 
which gives rise to bound Ai3-strings. Starting with a randomized high-energy field distribution, we found that an 
interconnected string network is indeed formed, consisting of A and B strings, as well as their (1,1) bound states, 
joined together at y-type junctions. No higher (p, q)-strings were observed in the simulations. 

Throughout each simulation, the network is dominated by a single ("infinite") net, with occasional small nets and 
a fair number of disconnected closed loops being formed in the course of the evolution. The characteristic length scale 
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FIG. 13: Fraction of total string length in bound strings, for flat, radiation and matter regimes, for simulations with n — 0.95. 
The percentage is still fairly low, though somewhat higher than for k — 0.90 (Fig. [5]) 



of the network approaches the scaling regime where it grows proportionally to time, ^(r) = ar with a ^ 0.15, in both 
radiation and matter eras, as well as in flat spacetime. Other characteristic length scales, such as the length of AB 
segments, the typical correlation length of A and B strings or the typical distance between Y junctions also scale. 

A surprising feature of our simulation is that bound AB-strings constitute only a small fraction of the total string 
length (~ 2%). Also, the average length of AB-segments is much shorter than the length of A- or _B-string segments. 
This is in contrast with analytic models (TRITsI predicting all lengths to be fairly equal. From movies of the simulations 
one can see that Ai3-segments do not always form when an A-string meets a _B-string (even for relatively high bounding 
energies); on the contrary, the formation of bound segments is rather infrequent. Even though new Ai?-segments are 
constantly formed, their lifetime is relatively short, usually less than the Hubble time, as the segments are "unzipped" 
by the free A and B ends pull them in different directions. 

Even though, with our general initial condition configuration, most of the typical distances in the network showed 
a scaling behaviour, the average comoving lengths of A and B segments did not. These lengths remained much larger 
than the typical correlation length throughout the simulation. This evolution regime cannot continue indefinitely, 
since the correlation length ^ would eventually catch up with the segment lengths and Ig. In an attempt to go 
beyond this transient regime and reach true scaling, we introduced an extra period of strong damping at early times 
(which has the effect of decreasing Ia/C) ^-nd pushed the parameters to extend the dynamic range of the simulation. 
We found that the system does seem to approach a regime where all characteristic lengths scale linearly with time. If 
this were the true scaling regime, it would provide us, amongother things, with the means to calculate CMB power 
spectra predictions from the field theoretical model as in f3l|, WK . IH, |4^ . But due to a relatively short duration of the 
linear evolution and to likely presence of spurious damping, our simulations cannot be relied upon for a quantitatively 
accurate description of scaling. 

Disconnected closed loops constitute about 10% of the total string length. Some of these loops are very large and 
will very likely reconnect into the main network again. If one factors out disconnected loops of length larger than 
a few times the correlation length, the total string length in loops drops to below 2%, in agreement with [4(|. We 
examined the length distribution of loops in the network and found that, even though out parameter choice is not 
expected to resolve accurately string dynamics, these distributions seem to scale. 

The network properties in our simulations are closer to those of superstring networks than they were in earlier 
simulations that used Za-strings or non-linear sigma-models. However, there is still an important differences. On the 
one hand, our A and B strings have the same tension, as opposed to F~ and D— strings. On the other, collision of 
same-type strings in our model always result in reconnection (unless the string segments are moving extremely fast 
[sot), while in the case of superstrings the reconnection probability is p < 1 and can even be small 0, [1, [5l|. This 
feature can be accounted for in Nambu-Goto and analytic models. 

The efficiency of various energy loss mechanisms by the string network remains a topic for future research. Energy 
loss to loop production appears to be substantial, considering that the length in loops at any time is a few percent of 
the total and that the loops do not stay around for long and rapidly decay. Another important energy loss mechanism 
in field theory simulations is direct particle emission from strings [i^ H^l- In fact, the analysis in [l^ shows that 
emission of particles and of tiny loops which immediately decay into particles is the dominant energy loss mechanism 
for a single C/(l) string network, so it probably dominates in our simulations as well. It has been argued in p^.l45ll46j 
that this effect is spurious and is due partly to insufficient resolution of the simulations and partly to excessive 
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amount of noise in the initial conditions. This issue is not completely settled, since some of the recent Nambu-Goto 
simulations [43j and analytic treatments ^53i] indicate continuous production of microscopic loops throughout the 
network evolution. 

In summary, what have we learnt from our simulations? We have demonstrated that an interconnected network 
of strings can indeed form at a symmetry breaking phase transition. This network shows no tendency to freeze to a 
static configuration. On the contrary, it appears to approach scaling, with all characteristic lengths growing linearly 
with time. Qualitatively, our results indicate that bound strings constitute only a small part of the total string length 
and that the A and B string segments are rather wiggly, having lengths significantly greater than the correlation 
length ^. The latter property leads to relatively frequent self-intersections and allows the network to lose a substantial 
fraction of its energy in the form of closed loops. Our simulations also indicate that the true scaling evolution may be 
preceded by a transient regime in which the comoving lengths of A and _B-segments remain nearly constant in time. 

Some of the shortcomings of our approach can be overcome in Nambu-Goto- type simulations (e.g., of the kind 
developed in [s^l) or in analytic models (along the hues of [13, [III). Either of these approaches, however, requires 
some microphysical input. For example, one needs to know under what conditions a bound string is formed in a string 
collision, what fraction of the binding energy of the newly formed string is radiated away, etc. The advantage of a 
direct field theory simulation is that it accurately represents the microphysics. A combination of all three approaches 
will probably be needed to reach a full understanding of network evolution. 
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